Unruh-DeWitt detector event rate for 
trajectories with time-dependent acceleration 



Luis C. Barbado^ and Matt Visser^ 

^ Instituto de Astrofisica de Andalucia (CSIC), Glorieta de la Astronomia, 18008 Granada, 
Spain 

^ School of Mathematics, Statistics, and Operations Research, 

Victoria University of Wellington, PO Box 600, Wellington 6I4O, New Zealand 

E-mail: luiscb@iaa.es; matt . visserOmsor . vuw. ac .nz 
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derivatives of the acceleration (jerk, snap, and higher derivatives). At the lowest order, 
we recover a Planckian spectrum with temperature proportional to the acceleration of 
the detector at each instant of the trajectory. Higher orders in the expansion involve 
powers of the derivatives of the acceleration, with well-behaved spectral coefficients 
with different shapes. Finally, we illustrate this analysis in the case of an initially 
inertial trajectory that acquires a given constant acceleration in a finite time. 
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1 Introduction 

The particle content of the state of a quantum field is an observer-dependent notion. 
The best known and most extreme example of this dependence is the case of Unruh- 
DeWitt detectors following constant acceleration trajectories in Minkowski spacetime, 
where one or several quantum fields are set to be in their natural vacuum states. These 
accelerated detectors, when coupled to the quantum fields, detect a particle spectrum 
with thermal Planckian shape and temperature proportional to their acceleration. This 
is the well-known Unruh effect, first proposed by Unruh in [1]. 

To probe the quantum state of a field, one has to use some form of particle detector. 
Probably the best known idealized model of a particle detector is the so-called Unruh- 
DeWitt detector [1, 2]. Formally, it consists of a point-like particle with at least two 
internal energy levels, coupled to a scalar field by means of a simple monopole (scalar) 
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interaction. One can interpret the probability of the detector increasing its internal 
energy by an amount hw, due to the interaction with the field, as the probability 
of the presence of a particle of this same energy hu in the quantum field. But the 
behaviour of a particle detector not only depends on the state of the quantum field, 
but also on the specific trajectory in spacetime followed by the detector. Thus, this 
notion of particle will also depend on the trajectory being followed by the detector. 
For example, if we make the detector follow an eternal constant acceleration trajectory, 
and set the quantum state for the field to be the vacuum state (vacuum as seen by 
inertial observers), we recover the thermal spectrum characteristic of the Unruh effect. 

For trajectories with (eternal) constant acceleration, and in general for (eternal) 
stationary trajectories, it is clear that the probability of excitation of the detector per 
unit of time is constant along the trajectory. On the other hand, when we want to deal 
with more general trajectories, several non-trivial issues arise. The quantity we wish to 
define is the probability transition rate at some instant of the trajectory. That is, for 
an ensemble of identical detectors following the same trajectory we wish to calculate 
the proportion of them that becomes excited per unit time around some instant of the 
trajectory. 

For a specific implementation of this proposal, one could use for example a switch- 
ing function controlling the detector interaction, so that the detector is coupled to 
the field only during some limited period of time. Excitations of the detector could 
then (loosely) be interpreted as particles detected during this particular period of time. 
However, in general the introduction of a switching function leads to spurious transi- 
tions (even for inertial observers) due to the switching process itself [3]. Instead, we 
will follow a somewhat different procedure, first introduced in [4] . Using the interaction 
picture, we formally calculate the probability of transition of a detector coupled to the 
field (in a constant time-independent manner), but up to some proper time r. Then, we 
differentiate with respect to the instant r to calculate what we shall call the response 
function, which is the quantity we are looking for. This quantity is well-defined, and 
finite, (at least for a wide class of 'non-anomalous' trajectories). We will also provide 
a clear and detailed physical interpretation of this response function in terms of realis- 
tic multicomponent detectors, further developing the interpretation pointed out in [5] . 
Nonetheless, the sharp cut-off implicitly generated when evaluating the probability of 
excitation only up to some specific instant, has the potential to introduce new regular- 
ization problems. These problems have been already discussed in the literature [3-6]. 
We will make use of a formula presented in [5] that avoids these problems, and that is 
particularly useful for our purposes. 

A particularly interesting regime for an Unruh-DeWitt detector is that in which 
it follows a trajectory with a slowly- varying acceleration g{T). One expects that in 



-2- 



this regime the detector perceives thermal radiation with a slowly-varying tempera- 
ture, proportional to the acceleration at each instant of the trajectory with the usual 
Unruh formula T(r) = h\g{T)\/ {2'KckB) [1]. We may call it the adiabatic regime. If 
we neglected derivatives of the acceleration higher than the first, the regime would in 
physical units be defined as the one in which cg'{T)/g{Ty ^ 1. (This is a constraint 
on the jerk in terms of the instantaneous acceleration.) Setting c — ?■ 1 as usual, this 
constraint is commonly rephrased as g'{T)/g{TY <^ 1 [7-9]. As explained in [9] (in 
a different and more general context), this is a physically reasonable requirement: If 
we want the detector to perceive a thermal spectrum with temperature cx g{T), the 
relative change in the acceleration g'{T)/g{T) during the time needed to detect a par- 
ticle of the typical energy in the spectrum, 1/ g{T), must be small. In this paper, we 
will attain this adiabatic regime precisely as a zeroth order approximation of a general 
asymptotic expansion of the response function of the detector, which we shall call the 
adiabatic expansion. This is an expansion in powers of the derivatives of g{T). As we 
said, a thermal Planckian spectrum with temperature proportional to the acceleration 
is obtained at lowest order. The other terms will be corrections to that result due to 
non-adiabaticity, and will become of higher importance the further away we are from 
the adiabatic regime. 

We begin in section 2 with a review of the definitions of the physical quantities 
involved in the description of an Unruh-DeWitt detector and their physical mean- 
ing, followed by a brief discussion of the regularization problems. We also generically 
describe a one- dimensional trajectory with a time-dependent acceleration. Then, in 
section 3 we use the expressions introduced in the previous section to calculate the 
adiabatic expansion of the response function, and in section 4 we discuss its relation to 
the high-energy behaviour of the detector. In section 5, we use the adiabatic expansion 
with an example of a time-dependent accelerated trajectory and discuss the results. 
We conclude in section 6 with a brief summary. 

2 General expressions for the transition functions 

Throughout the article, we will use natural units h = c = 1. Let us consider a quantized 
Klein-Gordon real massless scalar field 0(x), defined in a flat 3 + 1 spacetime back- 
ground, interacting with a localized detector D, which carries its own clock measuring 
proper time r, and which sweeps out some trajectory x{t) through spacetime. The 
quantum field has an associated Hilbert space T-Lp, and the detector has its own Hilbert 
space Hd- The combined system is defined on the total Hilbert space "H := 'Hf®'Hd- 
The quantum field (in the absence of the detector) has Hamiltonian Hp, and the de- 
tector (when decoupled from the field) evolves according to the Hamiltonian i^^i. 
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2.1 Transition probability 



Apart from the free evolution, the key point is that one postulates an interaction 
Hamiltonian Hi, which couples the detector to the field linearly. That is 



Hi{t) := ami^r) (p [x{t)] 



(2.1) 



where a is the coupling constant, and m(r) is the monopole moment of the detector. 
For the time being, we do not consider a switching function in (2.1), and thus the 
detector is switched-on during its entire trajectory. 

By assumption, we start in the asymptotic past with the field in its Minkowski 
vacuum state \iIjfo)i and the detector in its ground state |-E_do) of energy -Edo; so that 



(r ^ -oo)) = \i)FQ) ® \Edq) 



(2.2) 



We will work in the so-called "interaction picture" for clarity. The quantum time 
evolution operator is then given by the time-ordered exponential 



f/(r) = ^ 



exp ( i / drHj^r') 



(2.3) 



If the detector clicks, by jumping to some well defined excited state \Ed), with Eo > 
Edo, then this forces a collapse 



liPir)) It/jp) ® \Ed) 
The amplitude for this to happen is 



e^exp i / dr'Hjir'] 



I^Fo) ® \Edo) 



(2.4) 



(2.5) 



and the probability that the detector is found in its state \Ed) after an internal time r 
has passed, given that you are not looking at the final state of the field, is 



(2.6) 



where the sum runs over a basis of the Hilbert space Tip of the field. 

This probability, by definition, is a non-negative quantity smaller or equal than 
unity. During its evolution in r, it can do nothing but oscillate. However, one would 
like to have a quantity capturing the physical intuition that, during periods in which an 
acceleration is present, the detector would be continuously encountering new particles, 
producing some sort of cumulative effect. 
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2.2 Cumulative macro-detector 



Let us consider identical Unruli-DeWitt detectors, all of them prepared in the same 
initial state (in the manner described above) and following identical trajectories. Phys- 
ically, one might construct this by spreading the N micro-detectors into a tiny volume, 
so that the largest distance between two detectors is still much smaller that the mini- 
mum wavelength one wishes to measure. In this way, one can still consider the whole 
ensemble as pointlike. This would be a physical model for a more realistic macro- 
detector. Now, assume in addition that each micro-detector has a coupling constant 
e := a/y/N. Given a specific r, one can always take N large enough so that one can 
approximate the evolution operator by its first-order term in a perturbative expansion 
in the small constant e <^ 1. In this case we have the evolution for each individual 
component of the macro-detector 



mr)) 



i + i dr'Hjir') 



llpFo) ® \Edo) 



(2.7) 



where / is the identity operator, so that the transition amplitude can be written as 



\iPfq) ® \Edo) 



= ie{En\m{0)\EDo) f dr'e-'(^--^-o)-'(^^|0 [x{t')] \^fo)- (2.8) 

J — oo 

Here we have used equation (2.1) and the fact that the detector, left to its own devices, 
evolves according to 

m(r) = e"'^^" m(0) e'-^°". (2.9) 
Then, taking into account that 



(2.10) 



where J^^ is the identity operator in the Hilbert space Hf, the transition probability 
can be written as 

VEAr) = \{ED\mmEno)\' f dr" f dr' e-^^^-^^o )(-"-') W(r",r'). (2.11) 

J — oo J —oo 

Here, VV(r", r') is the Wightman distribution written in terms of the proper time of 
the trajectory, 

W(r",r') := (V^^o|0(x(r")) 0(x(r'))|^i.o), (2.12) 
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which can be exphcitly expressed as [10] 



>V(r", t') = 5 2-, (2.13) 

47r2 {[t(r") - t(rO - i£] - [x(r") - x(r')] } 

where the hmit e — )■ 0^ is imphcitly taken after integrating the function. 
As a final step let us define the quantity NE^ij) '■= ^'PEni'^)- Then 

^f^^{r) = \{ED\m{0)\EDo)f f dr" f dr' e-'^En-En,)(r"-r') y^(^"^^'). (2.14) 

J —00 J —00 

This can be interpreted as the average number of micro-detectors within the macro- 
scopic detector that are excited at time r. This quantity is no longer conceptually 
restricted to have a value smaller than one and in fact does exhibit cumulative effects. 

In the situations analysed in this paper the trajectories are always inertial in the 
past but are extended up to arbitrarily long time intervals. The previous formula can 
always be used under the assumption that in the limit AT — t- 00 (where AT is the 
observation time) the number of finally excited micro-detectors A/^;^ (which in general 
will grow with AT, possibly even going also to infinity), divided by the total number 
of micro-detectors A^, is such that Meo/N ^ 1. 

2.3 Response function (and click rate) 

The quantity we have calculated so far is the average number of detectors being in an 
excited state of energy Ee, at some time r. But the quantity that we are really interested 
in is the one that gives information about the amount of particles that "appear to the 
eyes" of an observer following a specific trajectory at a certain time, i.e. independently 
of the details of the detector. In this sense, we should note that the average number 
is the product of two factors with entirely different origins. The first one depends only 
on the characteristics of the detector. The second is the quantity [10] 

J^{uj,r):= r dr" r dr' e-'"^""-^') >V(r", r'), (2.15) 

J— 00 J— 00 

with u := Ee — Eeq- This quantity is "detector independent", and retains the cumu- 
lative property of the average number. (This function evaluated at r — )■ 00 is called 
the response function in Birrell & Davies [10], but we will reserve this name for its time 
derivative, as done in other sources [3, 5].) To better understand the behaviour of the 
detector at a certain moment of time one can calculate instead the time derivative of 
this quantity: 

n{u;,T):=—{u,T). (2.16) 
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Explicitly, this can be easily calculated to be 







n{u,T) = 2 3ft 



ds e"''^* >V(r + s, r) 



(2.17) 



— oo 



This is what we shall call the response function of an Unruh-DeWitt detector. As we 
can see, the quantity only depends on the past history of the detector, and thus is 
completely causal [4]. 

In the case of eternally stationary trajectories, either inertial or constant acceler- 
ation, this can be further simplified. We have then that the Wightman distribution is 
invariant under time translations along the trajectory, so that one can define 



This allows us to obtain the following particularly simple expression for this case: 



This expression fits the usual one for the response function dealing with a stationary 
situation, as is the case of a constant acceleration trajectory [10]. For this reason, one 
could be tempted to call this function click-rate function, characterizing the numbers 
of excitations produced per unit of time. However, this interpretation is not always 
completely appropriate as one can check that this function can sometimes attain nega- 
tive values, indicating de-excitation of the detector, typically in moments of decreasing 
acceleration. Nevertheless, as we will show below, its behaviour provides important 
insights regarding the time development of the detection process, and therefore in 
the particle perception process (see also the discussion in [5]). One could attempt to 
build more complicated models of more realistic detectors by incorporating notions of 
irreversibility and detector latency, but such refinements do not seem crucial to the 
questions we will explore. 

2.4 Finite-time detectors and the regularized Wightman function 

In this paper, we are considering a multi-part macro-detector that is always switched 
on, and we measure average number of micro-detectors excited at some time r. If we 
instead wanted to consider a detector with a different switching time or process, we 
should have written 



W(s) := W{s, 0) = W(r + s,t) = W(r + Ar + s, r + Ar). 




Furthermore 



W*{s) = W*(r + s, r) = W(r, r + s) = W(r - s, r) = W{-s). 



(2.19) 




(2.20) 




(2.21) 
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instead of (2.1) for the interaction term in the Hamiltonian, where ^(t') is now a 
switching function, that takes positive values during the interaction and vanishes during 
the periods of no interaction. This function should be integrated together with the 
Wightman distribution when calculating the response function. 

As pointed out in [3], the Wightman distribution is strictly speaking well defined 
only when integrated with a bump function, that is, a smooth function of compact 
support. In fact, this criterion is reasonable, as a physically realizable detector is never 
switched-on during an infinite amount of time, and it is never switched on or off in a 
(strictly) sharp manner either. We will call finite-time detectors the detectors with a 
bump function as the switching function. 

One can see that calculating the response function by using (2.17) corresponds to 
choosing the switching function 



where G is the Heaviside step function, i.e. the detector is on from the infinite past 
until a proper time r. This function does not fulfil any of the two requirements of a 
bump function. Thus, results obtained using (2.17) together with (2.13) are not strictly 
speaking correct. For example, Schlicht [4] showed that this procedure leads to non- 
Lorentz invariant results, giving non-zero detection rates even for inertial observers. 

However, in [3] it is also proven that we can approach with arbitrary precision 
the regime with an infinite switching time, and a sharp switch-off process, by using 
a concrete family of bump functions as switching functions. Consider a finite-time 
detector that starts switched-off in the asymptotic past [^{t') = 0], then it is turned 
on smoothly during a period 6 up to some time tq, remains switched on [^(r') = 1] 
during a period AT up to some time r, and is finally turned off again smoothly during a 
period 6. One can calculate the quantity this process, and with it the response 

function by differentiating with respect to the switching-off time r. This last quantity 
can be approximated by [3] 

"^finite-time ('^j t) = 



In the limit where AT — t- oo and S ^ 0, equation (2.23) approaches the following value: 



^(r') = e(r-r') 



(2.22) 



1 

2^ 





(2.24) 
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This is identical to equation (2.17), but replacing the Wightman distribution W by 
what we shall call the regularized Wightman function, which is defined by [3] 

This regularized quantity is well-behaved in the limit where the two arguments coincide 
r" — i- t'. If we compare (2.23) with (2.24), we explicitly obtain 

I ^finite-time (W, T") -^(w,r)| = 

1 , f cos(ws) 1] 1 6 



271^ J ^oo l[t(r + s)-t(r)]'-[x(r + s)-x(r)]' } 27c^AT \AT^ 

<^J^ + 0(4p^- (2-26) 



27r2AT V^^^ 

This result means that, insofar as the detector is switched on during a sufficiently long 
period of time (AT — )■ oo), and the switching timescale is sufficiently short ((5 — )■ 0), 
one can reach arbitrary precision using equation (2.24) for calculating the response 
function of this finite-time detector. Our calculations with (2.24) will then have this 
precise physical interpretation. 

Also, note that the quantity (2.25) is not merely a distribution but also well-defined 
as a function. It needs no £ — ?■ 0"*" limit calculation, as the pole at r" = r' is 'dodged' by 
subtracting the function l/(r" — r')^. It is thereby analytically much more tractable, 
and thus much more useful for our purposes. Note also that it is manifestly Lorentz 
invariant. 

2.5 Variable acceleration trajectories 

We are particularly interested in calculating the response function of a one-dimensional 
trajectory with time-dependent acceleration. For this class of trajectories, we can 
explicitly write the integrated equations for the coordinates (see for example [11, 12]). 
If we consider for simplicity a trajectory along the x-axis, we have: 

t{T) = to+ dr'coshx(r'), x{t) = Xq + drsinhx(r), y = Vo, 

J rn J TO 



Z = Zq. 

(2.27) 

Here 

X{r) ■.= Xo+ f dr'g{r'). (2.28) 



TO 



The function g{T) is the instantaneous proper acceleration of the trajectory along the x- 
axis, that is, as measured in an inertial and comoving reference frame, (to, Xq, yo, Zq) are 
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the initial coordinates at r = tq. Finally the initial velocity at tq is (cosh xo, sinh xo, 0, 0). 
Trajectories described in this way are also used in [7], where they also study the re- 
sponse of Unruh-DeWitt detector with time-dependent acceleration. However, their 
choice of the "time variable" with respect to which the response function is defined 
differs from ours and, actually, leads to a non-causal dependence of the transition rate: 
The detector, after having stayed in an inertial trajectory during its whole past history, 
could eventually click just because "it is prepared" to start to accelerate some time 
after the click. Because of this non-causal aspect to their calculation, their results are 
not directly comparable to ours. 

Inserting (2.27) and (2.28) in equation (2.25), we obtain the expression 



W{T + s,T;g] 



1 



J^ds'cosh(^J^' ds"g{T + s")^ ^ - ds' sinh (^// ds"^(r + s" 



(2.29) 



Here we have written the functional dependence on the acceleration g explicitly in 
the argument of the regularized Wightman function. Specifically, observe that the 
notation Vr(— , ; ] is a completely standard (but somewhat rare) notation designed 
to emphasise that this object is a function of the first two arguments, but a functional 
of the third argument. Note that none of the initial conditions of the trajectory are 
present in this expression. This is completely natural, as they merely correspond to 
the choice of the 4-velocity and origin of an inertial reference frame, and thus play no 
role in calculating the response function, which is Poincare invariant. 



3 Adiabatic expansion of the response function 

Equations (2.24) and (2.29) formally are enough to completely calculate the response 
function for an arbitrary acceleration history g{T). But except for very simple accelera- 
tion histories, because of the integrals involved in this expressions, one will not be able 
to calculate explicitly even the regularized Wightman function in (2.29). However, as 
announced in the introduction, we can use now equations (2.24) and (2.29) to calculate 
the adiabatic expansion of the response function in (2.24). This expansion will provide 
an arbitrarily good approximation to the exact value of the response function in the 
limit in which the acceleration varies slowly enough, in a sense to be precisely defined 
below. Note that, in order to calculate the numerical value of TZ{uj, r), we have to know 
the whole past history of the detector, via the acceleration history g{T'). As we will 
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see, within the framework of an adiabatic expansion we can appropriately calculate the 
response function just using local properties of the trajectory at some time r. With 
this aim, let us assume that the acceleration is a real analytic function 5'(t'), at least 
for t' < T. In order to obtain an asymptotic series expansion of the response function 
7^, let us introduce an auxiliary parameter a, and define a new auxiliary function W 
in the following way: Substitute 



g[T) ^g[a[T - r) + r) , 
in the formal expression of the Wightman function (2.29) so that 

1 

1 



(3.1) 



W {t + s,T,a]g\ :-- 



Jq ds' cosh ( ds"g(T + as"] 



Jq ds' sinh ( ds"g(T + as"] 



(3.2) 



It is obvious that the regularized Wightman function W given in equation (2.29) is 
obtained from this auxiliary function by setting a = 1. 

We can find a formal series expansion for (3.2) in powers of a. The coefficients 
of this series expansion will depend on s and on increasingly higher derivatives of g 
evaluated at r. Explicitly 

W{t + s, r, a; g] = 

Wo (s, gir)) + W, (s, (^(r), (?'(r)) a + W^ (s, (?(r), (?'(r), /(r)) a^ + . . . , (3.3) 



where 



and 



Wo (s, g{T)) := lim W {t + s, r, a; g] 

a— s>0+ 



W^4s,g{T),g\T),...,g^^\T)) 



1 d''W{T + s,T,a;g] 



(3.4) 



(3.5) 



a-s>0+ 



Once we have performed the expansion in powers of a, we formally evaluate the ex- 
pression in a = 1 to recover the value of the Wightman function in (2.29), (this is 
straightforward from (3.2)), but now written as an asymptotic series expansion: 

W{t + s, r; g] ~ Wo (s, g{T)) + W, (s, g{T),g'{T)) + W^ (s, g{T) , g' {r) , g" {r)) + .... 

(3.6) 
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Lastly, we will (again formally) calculate the integral in (2.24) term by term using (3.6). 
It is clear that each term in the series is integrable: First, for each value of s, the object 
W{t + s, T, a; g] is a real analytic function of a, and thus the coefficients Wn are finite 
for each s. Second all the coefficients tend to zero for s — ?■ — oo faster than because, 
for a small enough in (3.3), taking just a finite number of terms in the series in a should 
provide an arbitrarily good approximation to W{t + s, r; g], which actually vanishes in 
that limit faster than as can be easily checked. We then have 

7^(u;, r) ~ 7^o (u, gir)) + 7^l (u, gir),g'iT)) + 7^2 (u, ^(r), ^'(r), /(r)) + . . . . (3.7) 

where 

nn{u,giT),g'{T),...,g^^^'\r)) := 2 T ds cos{u s) {s, g^r), g' (r), g^^'\T)) . 

J —oo 

(3.8) 

Expression (3.7) is what we shall call the adiabatic expansion of the response function. 

Let us remark that expression (3.7) is not merely a short time expansion of TZ 
around r and, although closely related (as we will see later), it is also conceptually 
and mathematically different from a high-energy expansion in inverse powers of cu. In 
particular, the first four coefficients of the adiabatic expansion are 

noiu,giT))=\g{T)\^-^^ 



271 e^-^'^ - 1 ■ 



^1 {^^9{r),g\T)) =g{T 
T^2 {uj,g{T),g'{T),g"{T)) =|^(r)| 
^3 (w, g{T),g'{T),g"{T),g"'{T)) =g{T) 



9\r) 
9{ry 



9[r) 



9[r) 



9"\r) g\r)g"{r) 



9i^Y 



(3.9) 



Here u := u/\g{T)\, and the fij are dimensionless functions of the dimensionless 
variable u. The first term is simply the thermal spectrum for a temperature T = 
\g{r)\/27rkB ■ That is, as we already said, we recover a thermal Planckian spectrum 
detection as a zeroth-order approximation in the adiabatic expansion. Note that TZo 
depends only on the acceleration, while TZi also depends on the jerk, and 7^.2 also de- 
pends on the snap, etcetera. Moreover, TZq actually depends only on the absolute value 



-12- 



of the acceleration. This is what one should expect, as the sign of the acceleration g{T) 
by itself has no physical meaning in Minkowski spacetime. On the other hand, the 
term TZi depends on the sign of the quotient of the jerk by the acceleration, g'{T)/g{r). 
But this sign does have a physical meaning, indicating an increasing or decreasing ac- 
celeration. That is the reason for the TZq term to be multiplied by |(7(t)|, while the TZi 
term is multiplied by g{T). This sign dependency goes alternatively on and off with odd 
and even terms in the expansion, respectively, in the expressions (3.9). Apart from this 
alternating factor, the quantities between brackets in (3.9) are sums of dimensionless 
functions fij multiplied by powers of the dimensionless quantities {5'^'-*(T)/(?(r)'+^}. 
The explicit expressions for the first three fijiou) functions are 

h,,{u) =-^[l + 2uj F{uj) + oj^ Giu)] ■ 

uj {Cj [(3 + 47r2) + (-3 + 271^) cosh(27rc2;)] - 37r sinh(27rw)} _ 



/2,l(c^) = - 



24sinh*(7rw) 



/2 2(w) = r {-1 + 3^2;^ (-3 + 471^) + vrtl; [(-S + w^ (l -37r2)) coth(7rc:;) 

12sinh (ttcI;) 

-97rw (-2 + vrw coth(7rw)) sinh^^(7rw)] } , (3.10) 
where we have defined the functions 



F{uj) := / du ; G{uj) := / du ^ (3.11) 



oo 



Explicit evaluation of these integrals leads to formulae in terms of the digamma function 

il^i^x) ■= r'(x)/r(x): 

F{uj) = [^'(1 - iuj) - + iuj)] ■ G{Co) = \ \^"{\ - iu) + + ic^)] . (3.12) 

In figure 1, we plot the values of the first three fi,j{Cj) functions, together with the 
(normalised) thermal spectrum TZq {u,g{T)) /\g{T)\. Note that each term contributing 
to the adiabatic series has a well-defined shape in terms of the quantity u. Changing 
the value of the derivatives of g{T) only increases or decreases the overall contribution 
of each term (but not its shape) in the total output. That is, we can consider the 
total response function as a superimposed sum of different spectra. The first one is 
the thermal spectrum. The others have different shapes, and their contribution is 
proportional to different powers of {fi'^'^(T)/(7(r)'+^}. This might seem trivial, because 
actually we have done an expansion in derivatives of 5'(t). However, it is worth noting 
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Figure 1. Values of TZq (a;, ^(r)) /\g{T)\ (the normalized thermal spectrum, solid line), /i,i(a)) 
(dashed line), /2,i(w) (dotted line) and /2,2('^) (dash-dot line), as functions of a). 



that because of the way we have set things up each coefficient fij{uj) is a well-behaved 
spectrum in dj, without poles, and vanishing for w — )■ oo (no ultraviolet catastrophe). 

Finally, we have to say that, in general, we cannot guarantee the mathematical 
convergence of the adiabatic series, as we have formally integrated the series expansion 
of the Wightman function without any direct knowledge concerning its convergence. 
We can only be sure that this procedure leads to an asymptotic expansion in the 
dimensionless quantities {g^''\T) / g^rY^^}. 



4 High-frequency behaviour 

All the terms TZni^] 9{t)i 9'i^)i • • • ; 9^"'K^))^ ^^(^ thus all the functions fij{uj), tend to 
zero for u — )■ oo, as can be straightforwardly seen from their integral expressions (3.8). 
Let us study in some detail the behaviour of their high-frequency falloff. The inte- 
grals (3.8) that define the adiabatic expansion (3.7) have the following asymptotic 
behaviour for w — )■ oo (see [3, 13]): 

7^„ {u, giT),g'{T), . . . , g^-\r)) ^ 2 f -- j Wi'"^''^ (O, ^?(r), ^?'(r), . . . , g^^^r)) , 

m=l ^ 

(4.1) 

where, using (3.5), the coefficients Wn^ of this expansion can be written as 



Wi^^ (0, giT),g'iT), . . . , g^-\T)) = 1 £^ (r + r, a; g] 



■ (4.2) 

s=0 
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Now, note from (3.2) that the function W{t + s, r, a; g] has the following scaling prop- 
erty: 

Q^WiT + Qs,T,a;g]=WiT + s,T,Qa;Qg], Q e R. (4.3) 

If we use this scaling property for Q = —1 and take into account the invariance of the 
functional W with respect to the change g — > —g, we are led to the following result: 

W^^O;g{T),g'iT),...,g^^\T))=0, for n + I = 2p + 1, p e fi. (4.4) 

That is, any "odd combination" of derivatives with respect to s and a will cancel when 
evaluated at s = and a — )■ O"*". Noting that in (4.1) all the derivatives with respect 
to s are odd, we conclude that all the terms in that asymptotic expansion cancel for 
n even. That is, the even terms in the series expansion (3.7) will decay to zero faster 
than any inverse power of u, when a; — t- oo. Also, the odd terms will contain only even 
inverse powers of u. 

Also from the scaling property (4.3) we can show that the coefficients Wn'^ vanish 
for I < n. In order to obtain this result, let us first note that the order in which we 
take the derivatives with respect to s and a, and hence the evaluations at a = and 
s = 0, in (4.2) can be interchanged because we are considering real analytic functions 
gij), and thus W{t + s^r^a^g] is also analytic. Second, because of (4.3), we can see 
that d^W / ds''\s=Q is a polynomial in a of at most degree /. Indeed, 



d'-W (r + s,r, a; g] 




y + Qs, r, Q 


■,9] 
















ds^ 




s=0 


Q=l 






dW (r + Qs 


, ^, a; g] 




Ql 


' 1 


W{ 


T + s,T,Qa;Qg] 








Q=0 




0" 


Q=0 


s=l 



• (4.5) 



(For dimensional consistency in this calculation, consider a quantity T with dimensions 
of time, and then in the numerics set it to T = 1.) Hence the result that Wn^ vanish 
for I < n follows directly noting that (4.2) is proportional to the nth derivative of a 
polynomial of degree /. 

This means that the expansion of 7^„(a;; (^(r), (^'(r), . . . ,5f*^"^(r)) in inverse powers 
of uj will start from the term with uo~'^''^'^^\ The adiabatic expansion up to order n 
gives us automatically the high-energy behaviour up to order n + 1. This relation 
is natural, because higher energies should depend crucially only on what happens in 
smaller intervals of time. And the slower the evolution of the acceleration is, the better 
a short interval of time "represents" what happens in larger intervals. This conceptually 
connects the adiabatic expansion with the high energy expansion. In particular, up to 
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fourth order the high energy behaviour is exphcitly 



7^3 {u,g{r),g'{r),g"{r),g"'{r)) ~ i\g{r)g"'{r) + \g'{r)g"{r)^ u'^ + 0{uj)-\ 

(4.6) 



127r2 » 5-^ ^ ' 2' 



Summing, we see 

+ ^ - ^g{r)g"'{r) - lg'{r)g"{r)^ uj-' + 0{uj)-'. (4.7) 

The first term in the expansion fits with that found in [3, 5], while in [6] they aheady 
find terms up to order six, which also coincide with the terms in (4.7). 

We have to say one final word about the high-energy expansion. In [3] it is shown 
that, for a finite-time detector as described in subsection 2.4, the response function 
falls off faster than any power of u when u ^ oo. Thus, it seems that the high energy 
expansion we have found here is just an artifact of an infinite switching-on time, and 
a sharp switch-off. This is true if one considers the power expansion as the tendency 
of the response function for arbitrarily high energies. As it is already mentioned in [3], 
the inverse power behaviour will always be cured with a much faster falloff in this limit, 
if the detector is switched carefully. But one can still use the high energy expansion 
just as a tool for calculating the response function for a fixed u with arbitrarily good 
precision. Recalling considerations done in subsection 2.4, as far as one makes sure that 
the effects due to the switching processes are negligible with respect to the precision 
desired, the use of the series is perfectly justified, and does give the response function 
for a finite-time detector. 



5 Numerical results 

We will now consider a particular case for the acceleration function g{T), given by 

(5.1) 

The detector begins with no acceleration in the asymptotic past, then starts to increase 
its acceleration smoothly and, after a transient period of about At in proper time, it 
asymptotically tends to reach a final acceleration that we will normalize to be equal to 1 



g{r) 



1 + tanh 



At 
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(which, together with h = c = 1, determines a natural system of units for the problem). 
Note that this parameter At is closely related to the (inverse of the) parameter a used 
in order to obtain the adiabatic expansion (3.7). 

For this particular acceleration history the regularized Wightman function W can 
be found analytically with the help of symbolic manipulation software (such as Mathe- 
matica). On the other hand, the response function needs to be calculated numerically, 
something we achieve with very high precision, thus providing a very good profile with 
which we can compare the various expressions obtained in this paper. 

5.1 Numerical response functions: slowly switching on acceleration 

Let us describe some numerical results. First, consider the case At = 1000, so that the 
acceleration is very slowly varying (in comparison to the timescale set by its asymptotic 
value, of the order of 1). In figures 2 to 5, we plot some graphs of the response function 
TZ{u, r) for different values of u. We also plot what would be the response function if 
the detector were just detecting a thermal bath with temperature proportional to g{T), 
(the zeroth order approximation, as it is discussed in section 3), and the result when 
using the adiabatic approximation up to the third order. 

First of all, we see that, for all the energies considered, the fit with the adiabatic 
expansion up to the third order is virtually perfect. On the other hand, for low energies 
the match with the thermal detection seems to be perfect, while for higher energies 
the results start to differ. This might seem surprising when taking into account the 
results found in subsection 4: The adiabatic expansion up to order n coincides with 
the high energy expansion up to order n + 1, so that the higher the energies that we 
measure, the better the thermal approximation will be. In particular, the zeroth order 
adiabatic approximation must be exact for arbitrarily high energies. And this is true, 
simply because both the exact result and the approximation are actually zero. But 
it is true for the absolute difference between the exact value and the approximation. 
What figures 2 to 5 make evident is the relative difference, in comparison with the 
magnitude of the exact result (see the values in the vertical axis). And although the 
absolute difference between the numerical results and the thermal curve for arbitrarily 
low energies (figure 2) is greater than for energies w ~ 3 (figure 5), in a relative sense 
the difference is negligible for low energies, while for high energies it is much more 
important. Moreover, for energies still of the order of a; ~ 3, it is the contribution of 
the zeroth order term to the total spectrum which is actually negligible with respect 
to higher order contributions. In particular, the "bell shape" that appears in that 
figure corresponds, to a very good approximation, to the first term in the high energy 
expansion (4.7), proportional to g{T)g'{T). 
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Figure 2. Response function as a func- 
tion of T for u = 0, At = 1000. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line (fully covered by the solid one) 
the third order adiabatic approximation. 
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Figure 3. Response function as a func- 
tion of r for w = 1, At = 1000. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed (almost fully covered by the solid 
one) line the third order adiabatic approx- 
imation. 



1.4x10"' - 
1.2x10"' 

1.x 10"' 
8.x 10"" 
6. X 10"' 

4.x 10"' 



-1000 -500 



500 1000 1500 2000 



Figure 4. Response function as a func- 
tion of r for oj = 2, At = 1000. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 



Figure 5. Response function func- 
tion of r for uj = 3, At = 1000. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 



Finally, note also that, for r — ±oo, the three curves plotted in all the figures 
converge to each other. This is to be expected, since in this region the acceleration 
approaches a constant value, and thus the zeroth order adiabatic approximation must 
tend to the exact result. 
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5.2 Numerical response functions: rapidly switching on acceleration 
We will now plot the response function for the case At = 5 (figures 6 to 9). 




-0.01 - 
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T 

Figure 6. Response function as a func- 
tion of T for uj = 0, At = 5. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 
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Figure 7. Response function as a func- 
tion of r for UJ = 0.5, At = 5. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 




Figure 8. Response function func- 
tion of r for Lo = 1, At = 5. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 



Figure 9. Response function as a func- 
tion of r for Lo = 2, At = 5. Dots are 
the numerical solution, the solid line is the 
thermal spectrum approximation, and the 
dashed line the third order adiabatic ap- 
proximation. 



In this case, we are much further away from the adiabatic regime. In contrast 
to the previous case, even for arbitrarily low energies, deviations from the thermal 
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response function can be clearly seen. The adiabatic approximation (up to third order) 
starts to fail for lower energies, but it is still virtually exact for higher energies. Here 
differences with the third order approximation are bigger for lower energies even in a 
relative sense. This is again in agreement with the fact that adiabatic expansion gives 
"for the same price" the high energy expansion. Note in particular that in figure 6 the 
third order expansion seems to be a worse approximation than the zeroth order. This 
is what one must expect when staying far from the adiabatic regime, as we are dealing 
just with an asymptotic series, without (in general) any guarantee of convergence. 

6 Summary 

In this article, we have first calculated the transition probability for an Unruh-DeWitt 
detector in Minkowski spacetime, when interacting with a Klein-Gordon massless real 
scalar field in its Minkowski vacuum state. The detector is set to be in its ground state 
in the asymptotic past, and we calculate the transition probability when it follows a 
trajectory x(r') up to some time r. We then consider a cumulative multipart macro- 
detector, as an ensemble of many Unruh-deWitt micro-detectors with suitably small 
coupling constant. This allows us to use first-order perturbation theory to calculate the 
average number of detectors that are excited at some time r. Ignoring the detector's 
internal details, and taking the derivative with respect to the time r, we calculate what 
we have called the response function TZ{u, r) at some time r and some energy u of the 
detected particle. We also justify why, and up to what precision, we can consider an 
idealized situation with an infinite switching time (from the remote past) and a sharp 
switch-off of the detector at r when calculating the response function. 

Using the general expression for accelerated trajectories in one dimension, we man- 
age to write the response function as a functional of the acceleration g{T'). Then, we 
expand in powers of a formal adiabaticity parameter, which is inserted by hand to 
control how fast the acceleration changes. After calculating the expansion in terms of 
this formal paramter, we remove it from each term in the expansion. The zeroth or- 
der term is the thermal spectrum with temperature proportional to the instantaneous 
acceleration. That is the well-known Unruh effect result. Higher-order terms involve 
increasingly higher-order derivatives of the acceleration function, (jerk, snap, etcetera). 

We see that all the individual terms in the adiabatic expansion tend to zero in the 
high-frequency limit. There is no ultraviolet catastrophe for any of these terms. We 
also see that the even-order terms decay to zero much faster than any inverse power 
of u, while odd terms of order n in the adiabatic expansion are of order u~^'^~^^\ In 
other words, the adiabatic expansion to order n automatically gives the high-frequency 
expansion to order n + 1. 
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Finally, we use the adiabatic expansion with a simple but representative example 
of switch-on function for the acceleration. The numerical results obtained reflect most 
of the aspects we have discussed in the article. This example shows that the adiabatic 
expansion is a useful tool for generalizing the study of the Unruh effect to trajectories 
with variable acceleration. Further generalizations could include (among many other 
possibilities) more than one-dimensional trajectories, non-monopole coupling with the 
field to study possible anisotropics in the radiation, or the same problem for a field 
theory in a curved spacetime. 
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